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This paper discusses Skyrmions on the 3-sphere coupled to fermions. 
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1 Introduction 



The Skyrme model is a classical field theory which describes atomic nuclei 
and their low energy interactions It is a non-linear theory admitting 
topological solitons, so-called Skyrmions, labelled by a topological charge 
B which can be identified with the baryon number. When the theory is 
quantized the Skyrmions give a good description of nuclei, A resonance 
P] and even bound states of nuclei [HI IS [3]. A striking feature is that 
Skyrmions can be quantized as fermions [HllZl- When the theory is coupled to 
a fermion field two different ways of describing fermions appear in the same 
theory. Physically, the fermion field can be thought of as light quarks in the 
presence of atomic nuclei 8 . Hiller calculated the spectrum of the fermions 
in the background of a Skyrmion [Hj. However, the approximation for the 
Skyrmion configuration was very crude. The spectrum shows a spectral 
flow of the eigenvalues depending on the value of the coupling constant. In 
particular one mode crosses zero. This general behaviour has been verified 
by other authors ^1 ^2 ^] • To gain a better understanding we consider a 
generalization of the model. Since static field configurations in the original 
Skyrme model in flat space are topologically equivalent to maps — > S^, 
the model can be generalized to the base space being a sphere of radius L. 
The original model is recovered in the limit L ^ oo. It turns out that for 
L = 1 and baryon number B = 1 the Faddeev-Bogomolny bound is satisfied 
and the minimal energy solution can be worked out explicitly [El- In this 
case the map — > 5^ is an isometry. For larger L the Skyrmions are 
localized and tend to the usual flat space Skyrmions in the limit L — > oo. 

Monopoles and Skyrmions have many properties in common, and it is 
therefore instructive to recall the results for monopoles. There is an index 
theorem making use of the Bogomolny equations which states that for n 
monopoles there is a n dimensional vector space of fermionic zero modes jl4j . 
Furthermore, the monopole and its fermionic zero mode can be calculated 
explicitly for n = 1. We will examine how far we can extend the analogy 
between monopoles and Skyrmions. 

In this paper we study fermions coupled to Skyrmions on S^. The Dirac 
equation on S'^ is derived in Sect. |2 We also explain the Skyrmion- fermion 
interaction. In Sect. |21 we derive the ordinary differential equation which 
describes the solutions with total angular momentum G = for general 
shape functions. In Sect. 14.11 we give a brief discussion of ordinary differen- 
tial equations and their regular singular points. In Sect. 14.21 the equation 
of Sect. 131 is solved for a constant Skyrme configuration, which is equivalent 
to a mass term for the fermions. We also discuss the boundary conditions 
which are necessary to derive the spectrum. In Sect. 14.31 we derive the 
eigenfunctions and the spectrum for the hedgehog configuration with the 
special shape function f{fi) = fi. In Sect. El the results of Sect. 14.21 and 
14.31 are verified numerically. The spectrum for more general shape functions 
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is also evaluated numerically and compared with results in the literature. 
We find an interesting relationship between the topological charge of the 
Skyrmion and certain singular points. We end with a conclusion. 



2 The Dirac Equation on 

In this section we derive the Dirac equation on a 3-sphere of radius L = 1. 
First we define an atlas of two charts Xi and Xi which are the stereographic 
projections from the north pole N and the south pole S, respectively. Next 
we define vierbeins on both charts which uniquely determine a connection 
one-form using the metric compatibility condition and the torsion-free con- 
dition. The connection one-form induces the spin connection so that we can 
finally write down the Dirac equation for each chart. ^ At the end of this sec- 
tion we discuss the interaction term which couples fermions and Skyrmions. 

In order to construct the stereographic projection let be embedded in 
Euclidean with coordinates {xi,X2,X3,w). The points of can now be 
labelled by the projection from the north pole onto the plane through 
the equator introducing the coordinates Xi. Alternatively, one can project 
from the south pole S introducing the coordinates Xj. These two charts 
are well defined everywhere apart from their projection point A^ and S, 
respectively. In terms of the Cartesian coordinates we obtain 

X, = -^ and X, = -^. (2.1) 

I — w 1 + w 

In the following, we derive the equations of motion for both charts. As 
we shall see in Sect. 14. 2| this enables us to discuss the regularity conditions 
at the north and south pole in a simple way. We will also use the two charts 
to discuss a nontrivial symmetry of the energy spectrum in Sect. 15.21 

The induced metric on can be calculated from the Euclidean metric 
in R^. Taking into account that our base space is a Lorentzian space time 
M X S^, the metric can be written as 

9 = diag (l, , where = Xf + X| + A|. (2.2) 

The other coordinate system gives rise to the same metric replacing Xi by 
Xi. On the overlap the two charts are related via the transition functions 

X, = ^Xi. (2.3) 
Therefore, the coordinate vectors and 9^. are related by 
d^^ = {R^6^j - 2XiXj) dx, . (2.4) 

good discussion of frame fields and spin-connections can be found in 15 , and also 

in 
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Note that this transformation has the determinant det(-|^) = —R^. We 
now choose non-coordinate bases = Ca^dxu, and Cq = ead<? . It is 



convenient to choose diagonal vierbeins 6^'* and Ca such that 

PnO — 1 P — 1 

and (2.5) 
pi — _1±^ z i — l±M 

where all other components vanish. The minus sign occurs because Xi and 
Xi are related by inversion which changes the orientation. The vierbeins 
satisfy 

ea^ep^gf.i, = r]ag, (2.6) 

where rjap = diag(l, —1, —1, —1) is the flat Minkowski metric. The same 
relations are true for the vierbeins Ca- The non-coordinate basis vectors 
are related by an S'0(3)-transformation 

^i = ^{'2XiXj-Sij)ej. (2.7) 

This is in fact a 7r-rotation about the vector With this choice of vierbeins 
the connection one-form can be calculated from the metric compatibility 
condition uja/s = —t^/Sa and the torsion-free condition 

dr + o;"^ A 9^^ = 0, (2.8) 

where = e^judX'^ is the dual basis of e^. After a short calculation we 
obtain 

w"'^ = — ^ (x'^dXl^ - X^dX") , (2.9) 
1 -|- it \ / 

and the same relation holds for expressed in Xj-coordinates. The spin 
connection is now given by 

0^dX'^ = -ia;"/^S,;3, (2-10) 

where = | [7^,7/3], and the 7-matrices satisfy the standard anticom- 
mutation relations {70,7/3} = ^VafS- The Lagrangian for fermions in curved 
space time is given by 

Acrmion = ? (i^^e^^ (5, + 0«)) (2.11) 

With our choice of coordinates and vierbeins we obtain 

>Cfermion = ? {Xi, t) (^ij^'dt - H " ^i) ) (^i, , (2-12) 
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using the coordinates Xi and 




Aermion = ^{Xi,t) i-f'^dt + if ^^5^, " ^0 ' (2-13) 



for the coordinates Xi. Equation H2.13() can be transformed into (|2.12p by 
using Xi = -^Xi and the induced spinor transformation p of the fermion 
field iplXijt) which we calculate in the following. The spinor transformation 
p is defined as 

i>{Xi) =pij{Xi). (2.14) 

The transformation /9 is a lifting of the local Lorentz transformation T of 
the vierbeins in (|2.7j) . Note that the transformation Tij = {2XiXj — 6ij) can 
be rewritten as 

T{a)ij = cos a 5ij + sin a eij^X^ + XiXj{l — cos a) (2-15) 

for a = TT. Here Xj = Since we are only concerned with the spatial part 
we work with the Euclidean 6ij rather than with T]ij. Infinitesimally, the 
rotation is given by T{a)ij = 6ij + aeijkXk- This can be lifted to a spinor 
transformation by p{a) = exp{—^aeijkXkT,ij). For a = vr we obtain 

P = ilol5ljXj. (2.16) 

Note that there is a sign choice. As the spinor transformation p is an element 
of the double cover of the rotation group, the transformation —p is also a 
lifting of the rotation Tij. 

In the following we discuss fermions coupled to Skyrmions. The full 
Lagrangian C of this model is the sum of the fermion Lagrangian /2fermioru 
the Skyrmion Lagrangian £skyrmion and the interaction Lagrangian £int.- 
In this paper we consider fermions in the background of a static Skyrme 
configuration and neglect the backreaction. Therefore, we do not discuss 
the Skyrmion Lagrangian any further, and the interested reader is referred 
to ^7]. The interaction Lagrangian C\a\.. given by 

Ant. = -S'V' (<7 + i75T • TT) V', (2.17) 

where U = a + it ■ ir \s a, convenient parameterization of the Skyrme field. 
Equation (|2.17|) is the standard parity-invariant manner in which fermions 
couple to a Skyrme field and is also called the Yukawa interaction term. 
This Lagrangian was first introduced by Gell-Mann and Levy (without a 
fourth order Skyrme term), |18j . Also see ^2 fo'^ a discussion of various 
different chiral models. Recently, the Lagrangian (|2.17|) has been discussed 
by Balachandran et al. 8 . 

Note that tp \s a. spin and an isospin spinor. In our notation the Pauli 
matrices r act on the isospin indices, whereas the 7-matrices act on the spin 
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indices. It is convenient to split the Dirac spinor into two 2x2 spin isospin 
matrices "01 and ip2 such that 




(2.18) 



In this notation the 7-matrices can be written in terms of quaternionic 2x2 
matrices using i times the spin Pauh matrices a and the two dimensional 
identity matrix I2 as a basis for the quaternions. Since any complex 2x2 
matrix can also be expressed in terms of quaternions with complex coeffi- 
cients it is convenient to choose the same basis for the quaternions as before. 
Then the spin isospin matrices ipi can be written as = (^0^ ^2 + i'l'k (^k- 
With this notation the spin operators act on tp by left-multiplication, cjfc'f/'i, 
whereas the isospin matrices act on tjj by right multiplication 

Tki^i = ipicrl, (2.19) 
= -tpicr2crkcr2, (2.20) 

where the transpose in ()2.19|) guarantees the group structure, and 1)2.20(1 
gives a useful trick to evaluate af. 

As mentioned above, we make the approximation that the Skyrmion is a 
fixed static background field, i.e. we ignore the backreaction of the fermion 
fields. Furthermore, we only consider spherically symmetric Skyrmions. For 
-6 = 1, Skyrmions are believed to be spherically symmetric, but for B > I 
this is no longer true. However, we expect that some generic features of the 
Skyrmion- fermion interaction are not sensitive to the precise shape of the 
Skyrmion. For more details see Sect. 15.21 Spherically symmetric Skyrme 
fields are best written in polar coordinates 

U = exp (ififi) • r) , (2.21) 

where /(/i) is the "radial" shape function and is the normal vector in /i- 
direction, see ()3.5|) . The shape function /(^) can be calculated by solving the 
Euler-Lagrange equation of the Skyrme Lagrangian numerically. However, 
instead of using a numerical solution it is more convenient in this context 
to approximate /(/^) by an explicit function. Manton showed in jl9j that a 
good ansatz for L not too large is the conformal ansatz given by^ 

fill) = 2 arctan {k tan . (2.22) 

Here A; is a constant which has to be determined by minimizing the energy of 
the Skyrmion. This ansatz describes a Skyrme configuration which is either 
localized around one of the poles or uniformly distributed. The value k = \ 

■^For L — > 00 the ansatz tends to /(r) — 2 arctan Qualitatively, this is still a good 
ansatz. 
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corresponds to the 0(4)-symmetric ansatz of the previous section, whereas 
for A; ^ 1, or A; ^ 1, the Skyrme field is locahzed at the south pole, or the 
north pole, respectively. Manton also showed that for L < \f2 the stable 
solution is = 1. 

Finally, we can write down the Dirac equation for fermions coupled to a 
spherically symmetric background Skyrmion. It is given by 

[il^'dt - ii (^^^^. - - 3^^''^ ^(^-*) = 0' (2-23) 
using the coordinates Xi and 

f i7°5t + iy {\±^Q^-x\ -guA 4>{Xr,t) = 0, (2.24) 



using the coordinates X^. Here the Skyrmion gives rise to 

U^' = exp (^if{KR)Hr~^ , (2.25) 
and the same equation holds in the coordinate system X^. 

3 Solving the Dirac Equation for G = 

In order to solve the Dirac equation we introduce spherical polar coordinates 
for equation 1)2. 23(1 and ()2.24() . In terms of the coordinates x = (xi, X2,X3, w) 
in W^, the polar coordinates for the 3-sphere of radius 1 can be written as 

X = (sin /i sin 6* cos (j), sin// sin 6 sincp, sin /i cos 9, cos fj.) . (3.1) 

The angles {fi,6,(j)) can now be expressed in terms of Xi as 

arctan -j^zj + ^ for ii^ < 1, 
H = { ^ for i?2 = 1, 

arctan -^I^Y for i?^ > 1, 



arctan 



^/xTTx| 



(j) = arctan — -. (3-2) 
Xi 

In terms of the Xi, the relations for 9 and (p are identical to the ones above, 
and for the fj, coordinate we obtain 



arctan for -R^ < 1, 



2 

arctan -M- + yr for > 1. 



for i?2 = 1, (3.3) 
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/ — sin (j) \ 


and e, 




cos (j) 






K / 



Equation (tOaj) gives rise to the following Dirac equation in polar coordi- 
nates: 

+ (e,.d, - ^ - + ^) - g W^) = 0, (3.4) 

where = V'(A') t) ^"^^ the unit vectors (e^, ee, e^) are given by 

sin 9 cos (/ 
sin 6 sin 
cos 6 

U"'^ can be written as U"^^ = cos f{fj.) + i^^e^-T sin f{fi). Similarly, equation 
((T^ give rise to 

(z7°5, + ^Y {e,.d, + + - ^) - 5 U^^) ^ = Q. (3.6) 

The spin connection term could be removed by setting ih = — ^ — $ and 

_L COS ^ 

Tp = respectively. However, it is easier to discuss the regularity of 

the solution in terms of the original fields. 

We adopt the standard representation of the 7-matrices where 7*^ is 
diagonal. This is convenient because the Dirac equation is invariant under 
parity. 

Now, the spinor transformation p in (|2.16|) can be evaluated explicitly: 

With this choice of 7-matrices, we can derive the differential equation for the 
stationary solutions of the Dirac equation by setting il){Xi,t) = e^^^tl^{Xi): 

E^=( 9cosfip)_ a.p + ^ge,.Tsinfi^,)\ 

\ a p- ige^ • r sm / (^) -g cos / (^) J 

where 

a • p = -z (e,, • (5^ + ^) - O^e, • ade - sTTTSilT^e^ • ad^ . (3.10) 



The components of the matrix in equation 1)3. 9(1 commute with the total 
angular momentum operator G =L -|- S -|- I where L is the orbital angular 
momentum operator, the spin operator is S = ^o" and the isospin operator 
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is I = ^T. In polar coordinates, the orbital angular momentum L is given 
by 

f d d 
Lx = I \sm(j)— + cot9cos(j)— 
■f .9 ^ . , d 

Ly = I { -COSCj)— + COtOSllKp — 
\ 06 0(j) 

L. = -.|-. (3.11) 

Equation (|3.9|) is also invariant under parity P with 

P il:{Xi) = 70 V'( - Xi) and P X, p-^ = -Xi. (3.12) 

It is worth discussing the parity operator P in more detail. P acts on the 
7-matrices and the partial derivatives as follows: 



P7oP-i=7o and P 7, p-i = -7^, 
PdoP~^ = do and P di P-^ = -di. 



(3.13) 



Using the spinor transformation p in (|2.16|) we can show that 

P i^{Xi) = --foi^{Xi). (3.14) 

Therefore, the parity operator P which is defined for the coordinates Xi 
differs from P by a minus sign.^ 

The = 0"*" positive parity spin isospin spinors are given by 

V.(/x,^,<A)= f (l-cos/^)^G(cos/.)a2 y (3_^5^ 
y z (1 - cos /i) 2 F(cos fi) (e^ ■a)a2 J 

This gives rise to the following system of first order differential equations: 

where F, G and / are functions of p. Note that the G^ = 0^ spinor 

^(^,^,0)= f ii^-^o,^^)^F{co,p){e,.a)a, \ 
y (1 — cos ^) 2 G(cos ^)cr2 J 

gives rise to the same equation with g replaced by —g. 



^The origin of this sign difference is that the charts do not form an orientation atlas. 
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In Xj-coordinates positive and negative parity are interchanged such that 
the 0^ is given by 



1 ~ 



^P{^.,B,ct>)=^ ^(l + cos/.)^F(cos/.)(e,.a)a2 \ ^^^^^^ 

\ (1 + cos /i) 2 (^(COS ;U)(T2 / 



This gives rise to 



d„yGj [ -(E-gcosf) ^J^^-gsinf )\G 
Again the 0^ spinor 



(3.19) 



^(/x,e,c/>)= ( (l + cos/^)^G(cos/.)a2 \ ^320) 

y i (1 + cos fi) 2 F(cos fi) (e^ ■ cr) a2 J 

gives rise to the same equation with g replaced by —g. 

Equation ()3.16|) can be transformed into a second order differential equa- 
tion for G{u) where u = cos/i. Note that since /x G [0,vr] it is consistent to 
set sin ^ = \/l — cos^ //. 



(/-, 9\ (a 1 9 /'sin/Vl — \ d o r, „ 2osin/ 
(1 - u^) - 4u + 1 + ^^-^ — — — + ^2 _ 2 _ 9 ^ ; 
^ ^ du^ y E + gcosf J du ^ ^ 

, 3g/sin/ jl + u) g^fsin^f \ 

H — T-= T7 , = — — — qf cos / G(n) = 0, (3.21) 

2(^ + ffcos/) E + gcosf '^■^ JJ^' 'V ; 

and similarly equation 1)3. 19() can be transformed into 
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2\ / . 9 /'sin /a/1 — u'^\ d 2 q '^gsinf 



dn2 \ E-gcosf J du ^ ^ Vl - 

3qf'sinf (1 — u) f svo? f „, A x / x 

+?7^^ St -2 - V ^ + 5/'cos/ G{u) = 0. 3.22 

2[E- g cos f) Vl - V? E-gcosf J 

It is not easy to transform equation (|3.21|) and (|3.22|) into each other. In 
fact, due to the spinor transformation ()3.8() G{u) and F{u) are related via 



^(^) = Vr3^^(^)- (3.23) 

Also note that F{u) is given by 

F{u) = -—^ {{-^7=j+9^i^^f] G{u) - V^'^G{u)) , 

E + gcosf W 2V1 - n2 / du J 

(3.24) 

and similarly F{u) can be written as 

F{u) = i ( ( ~ -gsinf] G{u) - Vl-u^^G{u)] . 

E-gcosf \\2Vr^ y V ; V ^ V 

(3.25) 
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4 Explicit Solutions of the Dirac equation 



In this section we solve the Dirac equation which has been derived in the 
previous section for various shape functions /(//). However, first we review 
some basic facts about ordinary differential equations and regular singular 
points, following jJOj. 



4.1 Differential Equations and their Singular Points 

Any linear second order differential equation can be brought into the stan- 
dard form 

y"{u) + P{u)y'{u) + Q{u)y{u) = 0. (4.1) 

If P{u) and Q{u) are analytic functions in the neighbourhood of u = mq then 
u = uq \s called a regular point. In this case there are locally two linearly 
independent solutions yi{u) and y2{u) such that any solution of (|4.1|) is 
a linear combination of yi{u) and y2{y). Moreover, y{u) is an analytic 
function, i.e. it can be written as 

oo 

y(u) = Y,an{u-uoT. (4.2) 

n=0 

If either P{uq) or Q{uq) is singular then u = uq \s called a singular point. 
In the special case that 

= (4.3) 

U — Uq 

Qiu) = r^^,, (4.4) 
[u - Mo) 

where p{u) and q{u) are analytic functions, i.e. 

p{u) =po+pi{u-uq) + ... and q{u) = + qi{u - uq) + . . . , (4.5) 

the point u = uq \s called a regular singular point. Equations which only 
have regular or regular singular points are known as Fuchsian differential 
equations. Equation ()4.1() can then be solved by the following ansatz: 

oo 

y{u) = {u-UQY^an{u-UQ)''. (4.6) 

ra=0 

The exponent p can be calculated by solving the so-called indicial equation 

p2 + (po _ i)p + go = 0, (4.7) 

which is derived by inserting the ansatz (|4.6|) into (|4.1j) . The two solutions 
(/Oi,P2) of the indicial equation ()4.7() are called the exponents of the regular 
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singular point uq. Our convention will be that pi > p2- If the difference 
of the exponents s = pi — p2 is not a positive integer or zero, then the 
ansatz (|4.6|) gives rise to two linearly independent solutions of (|4.1j) which 
are uniformly convergent in some neighbourhood of uq. Note, however, that 
the solution diverges at uq if its exponent p is negative. 

Let s = pi — p2 be a positive integer or zero. This situation is slightly 
more subtle because the two solutions can become linearly dependent, and 
this could lead to further logarithmic terms. The solution yi{u) with ex- 
ponent pi is always uniformly convergent in some neighbourhood of uq. 
Therefore, we reduce equation (|4.1j) to a differential equation of first order 
by setting y{u) = yi{u)z{u). This gives rise to the following equation for 
z{u): 

{u - uofz"{u) + (2{u - uof^^ + (n - uo)p{u)] z'{u) = 0. (4.8) 

V yiw J 

This equation can be solved by integration and separation of variables, and 
we obtain 



W — Uq 



u 

= A + B2 ^ exp {-pi{v - Uq) - Ip2{v - uof - . . .) dv, 

J {yi{v)) 

u 

= A + B3 J{v- uo)-P<'-^P^ 9iv)dv, (4.9) 

where g{u) is an analytic function in a suitable domain, and A and Bi are 
constants. A short calculation shows that 

-po - 2pi = -s-l. (4.10) 

Let 

00 

g{u) = l + J29n{u-uor. (4.11) 

n=l 

Then for s 7^ the most general solution is 

y{u) = Ayi{u) + B {gsyi{u) \og{u - uq) + miu)) , (4.12) 

where A and B are arbitrary constants, and yi{u) is the solution of the form 
(|4.6|) . The constant Qs and the function y2{u) are completely determined by 
the differential equation 1)4. 8|) . The function y2{u) is given by 



\ ^ n=l 



y^{u) = {U- — + hniu - Uq)" , (4.13) 
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the coefficients hn being constants. In tlie special case tliat gs = no loga- 
ritlims occur, i.e. the two solutions yi{u) and y2{u) are linearly independent. 
If s = the corresponding form of the solution is 

y{u) = Ayi{u) + B (^i{u) \og{u - uq) + (n - uqY^ ^ hn{u - uq)"^ . 

(4.14) 

It is also possible to discuss the singularity at infinity. The transfor- 
mation V = maps the point at infinity to the origin. If we apply this 
transformation to equation 1)4.11) we obtain a new differential equation in 
terms of w. If the new equation has a regular or a regular singular point at 
the origin v = then (|4.1() has a regular or regular singular point at infinity. 

If a Fuchsian differential equation has n finite regular singular points 
u = Ui with exponents p\ ^ and a regular singular point at infinity with 
exponents then the following equation holds: 

n 

p(°°)+pr^ + ^(p«+4^))=n-l. (4.15) 

i=l 

This is a useful consistency check for the exponents. 



4.2 Fermions coupled to a Constant Skyrme Field 

We are interested in the behaviour of fermions in the presence of a back- 
ground Skyrmion. The Skyrmion introduces a space-dependent complex 
mass for the fermions. Therefore, it is instructive to consider first a massive 
fermion on S'^ . This corresponds to /(//) = so that the fermion mass is g. 
This example can be solved explicitly in terms of hyper geometric functions. 
In the following sections we proceed to more complicated configurations. In 
Sect. 14. 3| we consider the spherically symmetric shape function /(/u) = /i 
which can also be solved explicitly, but in a less standard way. In Sect. 
15. 2| the conformal shape function ()2.22)) is discussed. In this case, the dif- 
ferential equation can no longer be solved explicitly. We present numerical 
algorithms in Sect. 15.11 

Setting /(/i) = 0, equation ()3.21|) can be written as 

((1 - ^ - + 1) ^ + - i) ^(") = 0' (4-16) 

where e = \J E'^ — 9^ ■ 

Equation ()4.16)1 has three regular singular points, one at u = 1, one at 
u = —\ and one at u = cxd with exponents p\ \ = (0, — |), p{ 2^ = (0, — ^) 

and p^^^ = (| + e, I — e) respectively. Note that relation (|4.15j) between the 
exponents is satisfied, i.e. the sum of the exponents is equal to 1. In this 



13 



example we will show that the solution is non-singular on S*^ if its exponents 
are p^~^^ = and p^^^ = 0. 

Equation (|4.16|) can be solved explicitly: 

G{u) = ci.F(f + e,|-e;|;i(l + tx)) 

+ C2(l + 'u)"^.F(l + e,l-e;i;i(l + 'u)) , (4.17) 

where J^{a,b;c; z) is a hypergeometric function, e.g. see [211 P- 556]. The 
hypergeometric function is non-singular for < z < 1 and is normalized at 
2 = 0, i.e. T{a,b; c;0) = 1. It may be singular at z = 1, corresponding to 
u = 1, which is the north pole. As G{u) has to be non-singular at the south 
pole, u = —1, we have to set C2 = 0. 

Since G{u) is not well defined at u = 1 we have to express G{u) in terms 
of G{u) and calculate G{1). Using the transformation (|,3.8j) from ip to ip we 
obtain 

G{u) = (v— — ) {-i(T ■ ia ■ e^F{u) 



1 

l+u 
1 - u 



F{u), (4.18) 



which is the analogue of equation (|.S.23|) . Setting f{p) = 0, in equation 
(|H.22j) we can derive a similar equation for G{u) which has exponents p[ = 

(0, — |), p^l ~ (0; and p^^^ = (|+e, |— e). The solution of this equation 
can again be expressed in terms of hypergeometric functions: 

G{u) = ci.F(| + e,|-e;|;i(l-7x)) 

+ C2(l-u)"^.F(l + e,l-e;i;i(l-n)) . (4.19) 

This is only non-singular at u = 1 if we set C2 = 0. Assuming that E — g ^ 0, 
F{u) can be calculated from equation (|3.25|) : 

Therefore, G{u) can be written as 

^(^) = E^g + ")^^(") - • (^-21) 

G{u) has to be non-singular at u = 1. Equation (|4.21l then implies that 
G{u) is also non-singular at u = 1. Therefore, the regular solution G{u) has 
the exponents p^~^^ = and p^^^ = 0. Note that this argument only relies 
on the local behaviour and explicit solutions are not needed. As we shall 
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see, in the generic case, the exponents at the singular points do not depend 
on the chosen shape function, so that our analysis is also valid for general 
shape functions. Problems can occur if additional singular points coalesce 
with the singular points at « = ±1 because this changes the exponents.'^ In 
the generic case, the regular solution can be calculated by imposing that the 
solution can be expanded as a power series at the north and at the south 
pole. Later, this will be the basis for our numerical algorithm. Now, we 
impose this condition in order to calculate the spectrum of equation (|4.1(i)) . 

Using formula (15.3.6) in we can express G{u) in terms of hyperge- 
ometric functions at the second singular point. ^ 

G{u) = ci.F(| + e,|-e;|;i(l + 'u)), 

^ { 2 J r(|-e)r(| + e)-^^ 2'2i^ ^iiJ- 



The last term has to vanish, otherwise F{u) would be singular at the north 

;i + e) orr(| 



pole. Therefore, either r(| + e) or r(| — e) has to have a pole. Since e > 0, 



we obtain e = N + ^, and the spectrum is therefore 



E 



:Y''5^TcivTff, where TV = 0,1,2,... (4.23) 



With this condition the eigenfunctions GAr(n) reduce to Jacobi Polynomials 
(HU, formula 15.4.6). 

The spectrum is displayed in figure ^ It is symmetric with respect 
to £" = 0. The importance of this result lies in the fact that for g = 
the energy En and the corresponding solutions Gn{u) are the same for all 
possible shape functions /(//). Therefore, these solutions can be used as 
starting values for the relaxation method, which will be described in Sect. 
El As the radius of the 3-sphere tends to infinity the spectrum tends to the 
usual spectrum of the Dirac equation, which is continuous for E'^ > ff^- Also 
note that if we assume that G{u) and its first and second derivative remain 
finite as g tends to infinity then the energy E has to be of order g in this 
limit, so zizg is the leading order asymptotic for E. Therefore, we displayed 
the lines E = ±g in figure ^ As we shall see in Sect. [5] the region E'^ < g^ 
also plays a special role for our numerical algorithm. 

Another way of deriving the spectrum (|4.23|) is to note that if the ex- 
ponents at n = ±1 both vanish then G{u) is an integral function, that is it 

*In the following section we will also discuss the case that further singular points are 
inside the interval ( — 1,1). 

^Equation 14.221 only holds if |arg(l ~ z)\ < tv. However, this inequality is satisfied 
because (1 — = i(l + cos/i) > 0. 
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Figure 1: The energy levels En as a function of the coupling constant g for 
/(/i) = 0. E = ±g is shown as dashed lines. 

has only one singular point at infinity. The only integral functions which do 
not possess an essential singular point at infinity are the polynomials, e.g. 
[2{)\ p. 106]. Therefore, one of the exponents at infinity is equal to —N for 
= 0, 1, . . . . Since the exponents at infinity are p^^'' = (e + |, — e + |) we 
obtain the same spectrum as before. Note however, that since the exponents 
p^^^ differ by an integer, we have to check in this approach that the solution 
is indeed a polynomial, i.e. gg in equation (|4.11|1 vanishes. This spectrum 
has also been calculated by other authors, e.g. [22l l23]. 

4.3 Fermions coupled to the Skyrme Field /(/i) = fi 

In this section we construct the solutions of equation (|3.21l) for /(/i) = p. 
Configurations with f{p) / are called "hedgehog configurations" because 
at each point the vector in isospin space is proportional to the coordinate 
vector. For f{p) = p the map from the base to the target is even an 
isometry and possesses 0(4)-symmetry. Setting f{p) = p in equation (|:-{.21j) 
we obtain 



(1 - u^) {E + gu) 



{{^u + l){E + gu)+g{l-u^ 



))^+ (4.24) 
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{{E + gu) {E^ - + 2g - - gu) + lg{l + u) - g^ (l - n^)) G = 0. 

This is a homogeneous differential equation of Puchsian type with four sin- 
gular points. They are located atu = ibl,M = — ^ and u = oo. Assuming 
E 7^ zizg the exponents are 

(2,0), 

(^1 ± - + 4^2 + 8g- 4c,2^ . (4.25) 

For E = g the singular points at u = —1 and ti = — coalesce, and therefore, 
the exponents at u = — 1 become p[ — (^,0) ^^'^ the exponent at infinity 
is Pi^'* = (1 =b ^-y/l + Ag). The relation for the exponents (|4.15|) is satisfied. 
For E = —g the singular points at n = 1 and u = coalesce. This changes 

the exponents at u = 1 to p^l = (—1, — §)• In this case, there is an essential 
singular point at infinity and the equation is no longer Fuchsian. 

Since we are interested in solutions that can be expanded as a power 
series around u = ±1, it is natural to look for polynomial solutions. The 
singular point at u = — -j might impose a constraint. In order to verify 
whether there is a logarithmic term we have to calculate the coefficient gs in 
equation ()4.11|1 . Here s = 2. After a standard calculation we obtain g2 = 0. 
Therefore, there are two independent power series at the point u = 
This imposes no further restriction, since we do not have to exclude any of 
these two solutions. This implies that G{u) is again an integral function 
and we can calculate the energy spectrum using the exponents at infinity. 
Again, the exponents p^^^ differ by an integer, and we have to verify that 
the solutions are indeed polynomials. It turns out that the only problem 
occurs for p''^^ = 0. 

In the following, we derive the recurrence relation for these polynomials. 
Yet, it is instructive to consider first the solution where G{u) is a constant 
Oo, i.e. p^2°^ = 0. Then we can solve equation (|4.25|) for E and obtain 
E = g — ^ and E = ^ — g. Setting G{u) = uq in equation (|4.24j) . we obtain 

{E + g-^) (^g {E - g + ^) {u + 1) + {E - g + ^) {E - g)yo = 0. (4.26) 

It is easy to see that this is not a solution for E = g — ^. G{u) = ao is 
only a solution for arbitrary g if E = | — 5. This solution is very interesting 
because its spectrum crosses E = for g = ^. 



Pl,2 
Pl,2 

(00) 
Pl,2 
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Now we derive the general solution for equation (|4.24|) . Since we are 
interested in the solutions for which the relevant exponent vanishes for u = 
— 1, we expand G{u) around the south pole n = — 1: 

oo 

G{u)=Y,an{u + lT. (4.27) 

n=0 

Assuming that E ^ g, we derive the following recurrence relations. The 
initial condition is 

ai = -l(^E^ - {g-lfyo, (4.28) 
and the recurrence relation is given by 

_ [{E - 3g) n2 + {:iE -2g)n-{E- g) {e^ - (g - I)^)) 

" {E-g) (2n + 3) (n + 1) 

g[n^-[E^-E-g^ + 2g + \)) 
+ {E-g){2n + 2>){n + l) ^^'^^^ 

Since we have an initial condition we expect to be able to reduce this second 
order recurrence relation to a recurrence of first order. Therefore, we make 
the ansatz 

n-l 

a„ = 

Substituting this ansatz into the recurrence relation ()4.29() we obtain a non- 
linear first order recurrence for the 6„. However, this can be solved by setting 



in = n (4-30) 



{2E-2g + l + 2{n+l)){E'^-E + 2g-g'^ + \-{n+ 1) 

^" ~ {2E-2g + l + 2n) (2n + 3) (n + 1) 

(4.31) 

By construction the series (|4.27|) has the right behaviour at the singular 
point n = — 1. To find the behaviour at n = 1 we have to discuss the radius 
of convergence -Rconv. of the series. The radius Rcom. can be calculated as 
the following limit: 



-Rconv. — lim 

n— >oo 



lim 

n^oo \0r, 



1 



= 2. (4.32) 

Therefore, the series converges for u £ [—1, 1). For u = 1, equation 1)4. 27() 
can be written as 

oo 

G(l) = ^a„2". (4.33) 

n=0 
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For large n, 6„ behaves like ^. Therefore, a„ is proportional to (^)" so that 
the terms in the sum ()4.33|) are of order 1 and do not vanish. Hence, the 
series is divergent, unless it terminates. But because of the special form of 
Un and bn it is easy to see when the series terminates. In fact, for n > we 
can now write an as the following product: 



{E + g 



n\ (2n + 1)!! 



Y[iE^-E + 2g 



9' + \ 



(^ + 1) 



(4.34) 

where (2n + 1)!! is the product of all odd numbers less than or equal to 
(2n + 1). All an vanish for E = ^ — g which is the constant solution that we 
have derived previously. If the term (E + g— ^) vanishes then the coefficient 
for an vanishes individually. However, if one of the terms in the product 
vanishes, say for n = Nq, then not only a^o vanishes but also all ajy for 
^ ^ -^0) ^-e. the series terminates. 




Figure 2: The energy levels En as functions of the coupling constant g for 
/(//) = 11. E = ±g is shown as dashed lines. 



Prom this condition the complete energy spectrum can be calculated. We 
obtain the special level Eq = ^ — g which crosses from positive to negative 
energies, and two series of energy levels E^ and E^ which are given by 



n? + 2n+ [g — Vf' for 



n 



1,2,.... 



(4.35) 
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The spectrum is displayed in figure 121 If we only consider the modes with 
energy the spectrum is symmetric under reflection at the axes (7 = 1 
and E = ^. However, the mode with energy Eq is only symmetric under 
combined reflections g ^ 2 — g and S — > 1 — E, which corresponds to a 
TT-rotation around the point {E = ^,g = 1). This spectrum agrees with 
our previous calculation where we used the fact that G{u) is an integral 
function. 

5 Numerical Results 

In this section we first describe numerical algorithms to calculate the spec- 
trum for more general shape functions /(/u). Then we discuss the spec- 
trum for the conformal ansatz (|2.22|) . We also discuss the simple ansatz 
f{n) = Bfj, which gives rise to Skyrmions with B > 1. Finally, we compare 
our results with the literature. 

5.1 Numerical Algorithm 

The energy spectrum can be calculated for various shape functions using a 
standard relaxation method The differential equation (|3.21|) is rewrit- 
ten as a system of first order differential equations. 

yi'{u) = y2{u), 

y2{u) = h{yi{u),y2{u),E). (5.1) 

The energy E is considered as an additional dependent variable whose 
derivative vanishes, y3(tt) = E, so that the complete system becomes 

yi{u) = y2{u), 

V2{u) = h{yi{u),y2{u),y3{u)) , 

ys'iu) = 0. (5.2) 

For simplicity, let us discuss the numerical algorithm for 0(4)-symmetric 
shape function f{fi) = /x. Equation 1)4.24(1 has a regular singular point 
at u = —1 which implies that at this point we cannot choose value and 
derivative of the solution independently. However, once the value of G{—1) 
is chosen we can calculate G'{u) = ■^G{u). In fact, inserting u = —1 into 
equation (|4.24|) we obtain 

G'(-l) = -i {E' -g^ + 3g-l) Gi-1). (5.3) 

Similarly we obtain for u = 1, 

G'{1) = \{E^-g'+g-l + G(l). (5.4) 
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The system (j5.2|) is discretized such that the boundary conditions at the 
singular points u = ±1 are satisfied. Starting from an initial configuration 
the relaxation algorithm calculates an improved configuration by linearizing 
the equation around the solution and minimizing an error function which 
indicates how well the discretized equations are satisfied. 

The secret of a successful relaxation method is a good guess for the 
starting configuration. Fortunately, for g = the analytic solution of 1)3.21(1 
is known for all possible shape functions /(/^). Therefore, we start with this 
configuration at g = and increase g hy a small amount 6g. The relaxed 
solution is a good initial configuration for 25g, and so on. 

As a test, the spectrum for f{fi) = is calculated using the relaxation 
method. This reproduces the analytic result. For = fj- the relaxation 
method converges well as long as — ^ is not inside the interval [—1,1]. How- 
ever, if the regular singular point — is inside the interval, the relaxation 
algorithm no longer converges and has to be modified. This is one of the 
main reasons for including the lines E = zizg in the figures. 

In this case a specially chosen shooting method proves more suitable. It 
deals with the additional singularities in [—1, 1] by simply avoiding them. As 
mentioned before, our differential equation can be defined on the complex 
plane. Due to the boundary condition the solution G{u) is an analytic 
function. The original solution is recovered by restricting G{u) to be a real 
function of a real variable ti G [— 1,1]. In order to avoid further singular 
points on the real line we can set u = exp(z(/>) where (p £ [0,tt]. Then 
equation ([3.21(1 induces a differential equation for 0, that is the equation is 
solved on a semi-circular contour. The resulting differential equation can be 
written as a system of 5 real first order differential equations, one of which 
corresponds to the energy E. Now, the spectrum can be obtained by a 
standard shooting method, The boundary conditions for the shooting 
method are that at the singular points n = ±1 the solution is non-singular 
and furthermore that the equation is real at n = — 1 . Since the equation has 
real coefficients, the fact that G{u) is real for u = —1 implies that G{u) is 
real for u G [—1, 1]. 

A successful shooting method needs a good initial guess. For = we 
can rely on our explicit solution and impose G(— 1) = 1 and 



where equation (|5.5(1 follows from equation ((4.22(1 and equation ((5.6(1 follows 
from the reflection relation for the gamma function (see (6.1.17) in ^Ij). As 
in the relaxation method the initial guess for the coupling constant g + Sg 
is the value of the solution for coupling constant g. 

It is also possible to reconstruct the solution G{u) for u G [—1, 1] using 




(5.5) 



(5.6) 
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the solution G{exp(i(j))) for (p £ [0,7r]. Cauchy's formula states that if f{z) 
is an analytic function and C is a contour enclosing the point a then 



/(")(«) = ^/^i^^, (5.7) 

c 

see for example [201 page 90]. In our case G{u) is a real function on the 
interval [—1,1], so 

oo 

G{u) = Y,aku\ (5.8) 

A:=0 

where the coefficients are real. It follows that G{u) can be written as 



and a similar formula can be obtained for its derivative G'{u). For E'^ > 
both the relaxation method and the shooting method give the same result. 

5.2 Generalization to more complicated Skyrmions 

With the conformal shape function (|2.22)) . equation (|3.2H1 gives rise to a 
rather lengthy differential equation with five regular singular points, see 
|25j for more details. Assuming that the singular points do not coalesce we 
obtain the following exponents: 

pSS = (O'-i)' 



pif^^ = (2,0), 



pS = (0,0), 

p[f = (I±Ve^), (5.10) 



where a = k'^ + 1 and 6 = fc^ — 1 and k > 1. The relation (|4.15|) is again 
satisfied. Moreover, the exponents at infinity are the same as for /(/i) = 
0. In the limit /c — > 1 the exponents at infinity change discontinuously 
(see (|4.25j) ). The exponents change if one or more singular points coalesce. 
However, in general, there are five singular points which all have different 
characteristics. 

The regular singular points at n = ±1 provide the boundary conditions 
for the regular solution, and are not very sensitive to the exact form of the 
shape function. They imply that the solution of this equation locally still 
is an analytic function of u. One might expect that the regular singular 
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point at infinity could be used to calculate the degree of the polynomial as 
in the previous two sections. However, for general fc, the solutions are not 
regular on the entire complex plane, so that the exponents at infinity do 
not contain any information about the spectrum. This is because there is 
a further singular point at u = ^ with exponents (0,0). Generically, the 
solution diverges logarithmically at this singular point such that n = ^ is 
a branch point, see equation (|4.14|) . The existence of polynomial solutions 
might be related to the spherical 0(4) symmetry of the problem, which is 
broken for k ^ 1. Note that the singular point at | = f^^y independent of 
E, and it can never be inside the interval (—1, 1) because (fc^ + 1)^ > (fc^ — 1)^ 
for k G (0, cxd). Finally, the singular point at ^^Zga depends on the energy 
E. Moreover, for certain values of E, k and g this point can be inside 
the interval (—1,1). This singular point could force a potential solution to 
diverge, such that it would no longer be well-defined. Therefore, we are 
interested in the local behaviour near u = -Sr— After a straightforward 
but long calculation using Maple, we obtain 52 = 0. Therefore, locally, both 
solutions are convergent power series, and no logarithmic divergencies occur. 




Figure 3: The spectrum using (|2.22|) for the shape function /(/x) = 
2 arctan (A; tan ^) where k = 5. E = ±g is shown as dashed lines. 

In the following we calculate spectra as functions of the coupling constant 
g for different values of k numerically. For = 1 we reproduce the spectrum 
for /()u) = /i of Sect. 14.31 Also note that k = 00 corresponds to the spectrum 
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for /(/i) = vr which agrees with the spectrum for /(^) = in Sect. 14.21 As 
k increases the spectrum moves towards the k = oo spectrum. At the same 
time, the energy Eq intersects E = for increasing values of the coupling 
constant g, such that for A: — > oo the intersection point also appears to move 
to oo. In figure 01 we display the spectrum for A; = 5. 

In the following we discuss the eigenfunction which crosses zero in more 
detail. We restrict out attention to the following Lorentz scalar Q which is 
defined as 

Q= f^^lj. (5.11) 



53 

A short calculation shows that ipip is a function of /x only. Therefore, we 
can integrate out the angular coordinates 6 and (j) and obtain: 

TV 

Q = I Qil^) d/i, (5.12) 


7r 

= 47r / ((1 - cos/i) + (1 +cos;u) G^) sinV d/x. (5.13) 





The density Q{fi) measures where the fermion field is localised. To obtain 
equation (|5.13|) we used equations (|3.15|) and (|3.2H|) . G{u) and G{u) can be 
calculated from 1)3. 21() and (|3.22|) . respectively. Recall u = cos/U. Expressing 
F{u) in terms of G{u) avoids having to worry about the denominator in 
1)3. 24[) which is rather difficult to handle numerically. Note however, that 
F{u) is non-singular for E + (7 cos / = 0. This follows from the behaviour 
at the corresponding regular singular point of equation l|3.21() . The function 
G{u) is related to G{u) via equations (|3.23j) and (|3.24|) . This implies 

(5(1) = — -. (5.14) 

In the following, we calculate Q(/i) for k = 1, 2 and 5 and compare it to 
the baryon density 

B{fi) = -f'if,)sm^f{ji). (5.15) 
vr 

which measures where the Skyrme field is localized.^ 

The expression for ^jJ can be calculated explicitly for k = 1. In this case 
G{u) = oq. Since we are interested in the zero crossing mode we have to set 
E = ^ — g and obtain 

Q{lJ.) = STragsin^^. (5.16) 

Surprisingly, this is independent of g. Also note that in this case Q{fi) is 
proportional to B{fj,). In the following, we normalize Q{fi) so that Q = 1. 

"Usually, -B(/i) is defined with a minus sign, but this is of no relevance here. 
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For A; > 1, the Skyrme configuration is localised near the north pole 
11 = (see figure 0]). In figure 0] (a) and (b) we show B{fi) for k = 2 and 
k = 5, respectively. In figure E](c) and (d) we show Q{fi) for several values 
of the coupling constant g. For g = 0, Q{p) can also be calculated explicitly 
and is independent of k. For k > 1 and g > 0, the fermion configuration is 
also localised near the north pole. For k = 2, the energy E vanishes when 
g ~ 1.795. For k = 5, the energy E vanishes when g ~ 3.573. In figure |^(e) 
and (f) we display Q{fi) for these values of g and compare it to the baryon 
density B{fi). There is quite a good agreement between these two functions. 
Therefore, the E = states can be interpreted as a fermion-Skyrmion bound 
states. 

The simplest method for obtaining a Skyrmion with baryon number B > 
1 is the ansatz f{fj.) = Bfi. For B > 1 this Skyrme configuration has very 
high energy. However, we are mainly interested in the interplay of topology 
and the spectrum of the Dirac equation which should be independent of the 
precise form of the Skyrmion. 

For B = 2 the spectrum is displayed in figure El It is symmetric with 
respect to the E-axis. Note that the symmetry of the spectrum can be 
derived from equation (|3.21|) and l|3.22j) . For even B the function cos / is 
an even function of u and sin/ is an odd function of u. Under g — > —g 
and u ^ —u equation (|3.2H) transforms into equation (|3.22j) . However, the 
equations are just expressed in different coordinate systems, and hence have 
the same energy eigenvalues. Therefore, we have shown that the spectrum 
is symmetric under g — > —g if B is even, provided that / is an odd function 
of u. 

For (7 = all the eigenmodes have the same value, independent of the 
background configuration, and parity. As g increases, the modes increasingly 
feel the infiuence of the background configuration, in particular its topology. 
It appears that the winding number B = 2 of the Skyrmion forces 2 modes 
to cross zero, one for g > and one for g < 0. In order to follow the 
spectrum for E"^ < g'^ the shooting method has to be used. 

This pattern can also be verified for B = 3, see figure El Now there 
are three modes which cross zero, two of positive parity {g > 0) and one of 
negative parity (^f < 0). It seems probable that the number zero-crossing 
modes of positive parity differs at most by one from the number of zero- 
crossing modes of negative parity. Further calculations for i? = 4 and 5 
confirm this pattern. For i? = 4 we find two zero-crossing modes of positive 
parity and two of negative parity, whereas for B = 5 there are three positive 
parity modes and two negative parity modes. 

This numerical evidence shows that the singular points of equation (|3.21() 
play an important role. In the following we will assume that /()u) is an 
analytic function, and in particular that it can be expanded in a Taylor 
series with respect to the variable u = cos /i, for n E (—1,1). It follows that 
cos / and sin / are also analytic functions of u. Therefore, the only singular 
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points of (|3.21|) inside the interval (—1, 1) occur for cos / = —E/g. If such 
a singular point is non-degenerate its exponents are and 2. Furthermore, 
it can be shown be shown by expanding / as a power series in u that the 
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Figure 4: In this figure we display B{ii) and Q{ii) as a function of /i for 
k = 2 and 5, and for various values of the coupling constant g. 
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Figure 5: The spectrum using /(/u) = 2/i. E = zizg is shown as dashed Unes. 



coefficient §2 always vanishes. In our previous examples f{fi) was always 
monotonic. Let f{fJ.) be monotonic with boundary conditions /(O) = 
and /(vr) = Btt. If < then cos/ = —E/g has B nondegenerate 
solutions for G (— 1, 1). In this case, these singularities of equation (|.S.21j) 
all have exponents (2, 0) and the coefficient g2 vanishes, so that the solution 
of (|3.2H) is regular at these points. We conjecture that the existence of B 
regular singular points inside the interval (—1, 1) indicates that there are B 
modes crossing from one side of the spectrum to the other. This conjecture 
is equivalent to the statement that for topological charge B there are B 
modes that cross zero.^ 

If we allow for more general Skyrme configurations than the hedgehog, 
such a statement becomes even more difficult to prove, because now we have 
to solve partial differential equations. However, both the topological number 
and the number of singular points are stable under small perturbations. 
Therefore, there seems to be a deep interplay between the topology of the 
configuration and the singular points of a differential operator. This can be 
compared to monopoles where we have an index theorem which connects the 

^The existence of such zero-crossing modes has also been discussed recently by Bal- 
achandran et al. in a more geometric language. The authors discovered a relation between 
these modes and Euclidean instantons. 
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Figure 6: The spectrum using /(/x) = 3/u. Here we only show the spectrum 
for \E\ > \g\, see text. E = zizg is shown as dashed hues. 



topology of the monopole configuration with the index of the Dirac operator. 

In the following we compare our results with the literature. Kahana 
and Ripka considered the problem of coupling a hedgehog configuration to 
fermions in flat space using an approximation for the Skyrme configuration 
with a size parameter. Varying this parameter they computed the spectral 
flow of the eigenvalues. They found a O"*" zero mode for a special value of 
the "size" of the hedgehog. Since the size is given by minimizing the energy 
density for the Skyrmion their result corresponds to varying the coupling 
constant g. 

In ^01 a linear approximation for the shape function, /(r) = ir{l — 
was used. Here x = g ■ r so that the Dirac equation for the fermions is 
independent of g. However, we defined the Skyrmion as a given background 
configuration which minimizes the energy. Therefore, X = g ■ i?sk. where 
i?Sk. is the size of the Skyrmion. i?sk. can be calculated by minimizing 
the energy (2^: i?sk. = 1-67. Kahana and Ripka showed by numerical 
calculation that there is a zero mode for X = 3.2 which corresponds to 
g = 1.9. In |E] an exponential function was used to approximate the shape 
function: /(r) = 7rexp(— -|^). In this case minimizing the energy results 
in i?sk. = 1-7 and the zero mode occurs for X = 2 which corresponds to 
g = 1.1. The fact that there is a mode which crosses zero seems to be 
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independent of the special ansatz for the shape function. This is consistent 
with our results. 

In [nj it has been argued that a Skyrmion in flat space is reasonably well 
approximated by the stereographic projection of a Skyrmion on a 3-sphere 
of optimal radius. The optimal radius for i? = 1 is L = 1 so we can use our 
analytic result for the mode Eq in Sect. I4.ISI For 5 = | this is a zero mode, 
and the value of g is in the same order of magnitude as in the literature 
mentioned above. 

6 Conclusion 

We have derived the Dirac equation for fermions on the 3-sphere which are 
chirally coupled to a Skyrmion. In this paper, we only considered spherically 
symmetric Skyrme configurations and fermionic states with the total angular 
momentum G = 0. We were particularly interested in the shape functions 
/(/i) = and /(/Lt) = fj. and the conformal ansatz 1)2. 22() . 

The shape function /(/i) = leads to a hyper geometric equation and 
can therefore be solved explicitly in terms of hyper geometric functions. To 
derive the spectrum we have to impose that the solutions are non-singular at 
the north and south pole using the correct chart. We derived the spectrum 



which agrees with the spectrum calculated by Sen |22] using Maurer Cartan 
forms. Moreover, the eigenfunctions G]sr{u) are given by Jacobi Polynomials. 

The shape function /(^u) = fi is more complicated. It leads to a Fuchsian 
equation with four regular singular points, two at the boundary u = ±1, one 
at infinity and one depending on E and g. However, the equation could still 
be solved in terms of polynomials. We were able to solve the recurrence 
relation for the coefficients of the polynomial and derived the corresponding 
spectrum: 




for 



iV = 0,1,2,... 




for 




The eigenfunctions have been evaluated explicitly: 




n 



Gn{u) 



y^Qfc (li+i)'', 



A;=0 



where 



{E + g-D {E + g 
k\ (2A: + 1)!! 



ll{E'-E + 2g-g' + l-i^ + lf). 
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We also considered more general shape functions numerically. For B = 1 
we discussed the conformal ansatz (|2.22() and calculated the spectrum. For 
k = 1 the spectrum agrees with the analytic solution. In particular there is 
a mode whose energy Eq crosses zero for a certain value of the coupling 
g. The graph of Eq as a function of g changes as the parameter k in the 
ansatz is varied. However, for finite A; > 0, this mode always stays present. 
According to our numerical evidence the eigenfunction for E = can be 
interpreted fermion-Skyrmion bound state. 

For larger B we considered the ansatz /(/x) = B/j,. This leads to B 
modes crossing zero. We conjectured a connection between zero modes and 
certain singular points in the equation which in turn are connected to the 
winding number of the configuration. 
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